Creep dynamics of elastic manifolds via exact transition pathways 
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We study the steady state of driven elastic strings in disordered media below the depinning 
threshold. In the low-temperature limit, for a fixed sample, the steady state is dominated by a 
single configuration, which we determine exactly from the transition pathways between metastable 
states. We obtain the dynamical phase diagram in this limit. At variance with a thermodynamic 
phase transition, the depinning transition is not associated with a divergent length scale of the steady 
state below threshold, but only of the transient dynamics. We discuss the distribution of barrier 
heights, and check the validity of the dynamic phase diagram at small but finite temperatures using 
Langevin simulations. The phase diagram continues to hold for broken statistical tilt symmetry. 
We point out the relevance of our results for experiments of creep motion in elastic interfaces. 



I. INTRODUCTION 



Disordered elastic systems arc ubiquitous in nature. 
They appear as interfaces such as magnetic [l|, Jl, [3, Q 
or ferroelectric domain walls, contact lines [31, frac- 
tures l8|,|9||, or as periodic structures, such as vortex lat- 
tices'^ charge density waves [l2| and Wigner crys- 
tals . They all share the competition between elastic 
forces that tend to order the system and the microscopic 
disorder that seeks to distort its structure. This compe- 
tition manifests itself in the static properties, leading to 
a roughness of the interfaces or a distortion of periodic 
order. 

In addition, disorder leads to pinning and thus also 
affects the dynamical properties. Indeed practically all 
the above systems may be displaced through an external 
force (magnetic or electric field for magnetic or ferroelec- 
tric domain walls, current for vortices, etc.). The mo- 
tion directly influences central observables of the system 
(e.g. magnetization for magnetic domain walls, voltage 
for vortices etc.). 

The static properties of disordered elastic systems are 
now well understood. It has for example been established 
that interfaces (on which we will focus in this paper) 
become rough in the presence of disorder. The roughness 
is characterized by an exponent Coq, which depends only 
on the universality class of the disorder, the dimension 
of the interface and the nature of the elastic forces. The 
dynamical properties are less well characterized. At zero 
temperature, disorder leads to the existence of a critical 
pinning force /c, the depinning threshold, below which 
the interface is immobile, and above which steady-state 
motion sets in. The dynamics at finite temperature is 



even more difficult to analyze than the zero-T behavior. 

It has been particularly fruitful to consider the depin- 
ning transition as a regular critical phenomenon, with the 
velocity playing the role of an order parameter 1J| . In 
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this framework, the depinning transition appears linked 
to the existence of a correlation length ^ which diverges 
at the depinning transition, and to the presence of criti- 
cal exponents, both for this length C ^ (/ ~ fc)""'^''^ and 
for the velocity w ~ (/ — fc)^- For zero temperature, the 
analogy of the dynamical depinning transition with equi- 
librium critical phenomena has been checked directly by 
numerical approaches and by analytical techniques such 
the functional renormalization group. 

How to extend the analogy with critical phenomena 
to finite temperatures has not been completely evident. 
Measurements of the thermal rounding of the depinning 
transition, or observations of scaling for transient dynam- 
ics indicate that these ideas carry over to finite tempera- 
ture [H, [H, • However a direct study of the motion of 
such systems is very difficult for forces below the depin- 
ning forces, at finite temperature. Indeed in that case, 
the motion takes place by thermal activation over bar- 
riers, leading to extremely long activation times. This 
renders numerical techniques such as the molecular dy- 
namics inefficient. 

At finite temperature, the dynamics of a disordered 
system is very difficult to simulate, because of the high 
excitation barriers and the presence of thermal noise. 
Standard dynamical algorithms become very inefficient 
because the system is frozen in a local minimum. Ad- 
vanced dynamical simulation methods, such as the BKL 
algorithm and its variants cannot always be applied 
because the systems are in fact equilibrated on some 
length scale, leading to a futility problem [l^, [23|: Not 
only does it take a long time to pass a barrier, but before 
doing so, the system will have undertaken a very large 
number of (futile) moves among local configurations. 

This problem can often be overcome in statics, where 
we need not follow the very slow dynamics because ob- 
servable averages are given by the Boltzmann weight 
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for each configuration. For this reason, many meth- 
ods (transfer matrix, optimization algorithms, advanced 
Monte Carlo methods) allow one to characterize the ther- 
modynamics of disordered elastic systems with an effort 
polynomial in the system size. 

Similarly, although analytical techniques such as the 
functional rcnormalization group have been used with 
success to investigate the motion at finite temperatures, 
they lead to complicated equations that have been solved 
so far only for the 'creep' regime of very small forces 
f 0. Moreover, these techniques rely on an expansion 
around four spatial dimensions and are thus not very 
well adapted to tackle quantitatively realistic one- or two- 
dimensional interfaces. 

In Ref [2l| we introduced a novel numerical method 
which allows to follow the motion of an interface at finite 
temperature, without running into the above-mentioned 
difficulties. This algorithm was used to demonstrate 
that the analogy between the dcpinning and critical phe- 
nomena is incomplete. In fact, the steady-state motion 
lacks the divergent length for / — *■ /c from below. In 
the present paper, we discuss the algorithm in detail. 
We use it to study the various properties of an inter- 
face close to dcpinning. In addition to the questions of 
the steady-state motion and the corresponding divergent 
length scales and the roughness of the lines, we study 
the distribution of activation barriers during the motion. 
This is crucial since disordered elastic systems are glasses, 
with a priori divergent barriers. We also compare the re- 
sults of our algorithm to molecular dynamics simulations, 
in order to check that taking the T — > limit before the 
thermodynamic limit does not introduce artifacts. 

The outline of this article is as follows: After a 
short review, in Section |TT1 of the statics and the zero- 
temperature dynamics of disordered elastic systems, we 
outline in Section IIIII the finite-temperature dynamical 
phase diagram obtained with our method. Section IIVI 
discusses the special properties of the low-temperature 
dynamics, which are used in the algorithm, and SectionlVl 
contains a detailed description of our numerical results. 
A general discussion, including the prospects for exper- 
iments in Section IVIl concludes this paper. Technical 
details of the algorithm and several mathematical proofs 
of key properties of the low-temperature dynamics are 
contained in Appendix 1X1 whereas Appendix IB] resumes 
several properties of long-range elastic systems. 



II. BASIC NOTIONS 

In order to render our paper self-contained, we review 
in the present section how the competition between dis- 
order and elasticity manifests itself in the statics and 
the zero temperature dynamics of elastic manifolds. We 
stress the difference between equilibrium properties and 
the non-equilibrium steady-state behavior. Furthermore, 
we introduce the three 'reference states', the corner- 
stones of our analysis of the dynamical phase diagram 




FIG. 1: (Color online) An elastic string h{x) in a random 
medium, (a) For random-bond disorder (RB), its interaction 
with the impurities is local, (b) For random-field disorder 
(RF), the energy _Edis[/i] depends on V{u,x) for all u < h. 




FIG. 2: (Color on line) Velocity-force characteristics of an 
elastic manifold at T = 0. The three self-affine reference 
states are: Equilibrium (/ = 0); Depinning (/ = fc); Fast 
flow (/ > /c). 

of Section Uni 

A. Elastic manifolds 

Wc consider a d-dimcnsional interface separating a 
d + 1-dimensional covering space into two regions. The 
interface is free of overhangs or loops, and may thus be 
described by a uni- valued displacement field h{x). Fig. [1] 
shows an elastic string (d = 1) in a two-dimensional ran- 
dom medium, the case that we concentrate on in this 
paper. 

The elastic energy -Eei[^] is minimal for the flat man- 
ifold h = const. For short-ranged elastic interactions, 
deviations from this configuration are often described by 
the harmonic energy 

Eci[h] {d.hf (1) 

(with an elastic coefficient c). This case applies to mag- 
netic and ferroelectric domain walls as well as to vortex 
lattices and charge density waves. 
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(at / = 0); the depinning (at / = /c), and the fast flow 
(at / ^ /c). In these states, the manifold is spatially 
self-afSne. This means that lengths x and displacements 
h above a small cut-off (set by the size of the correspond- 
ing Larkin domain [sij ) can be rescaled as x' = hx and 
h! = \fh into a new interface h'{x') which is statistically 



equivalent to h{x). ( is the state's characteristic rough- 
TABLE I: Zero-temperature critical exponents for an elas- ness exponent. Even at finite temperature, manifolds can 
tic string in the three self-afBne reference states (equilibrium, be analyzed in terms of these reference states, as we will 
depinning, fast flow). discuss later. 



Dipolar forces 0, 0, , contact lines for wetting [2^ , 
and also crack propagation [23| represent the class of 
manifolds with long-range elastic interactions. These sys- 
tems are described by different harmonic forms (see Ap- 
pendix |B]). 

A second essential interaction is provided by the cou- 
pling of the manifold to disorder. Two types of disorder 
correlations have been much discussed in the literature: 
Random-bond (RB) disorder corresponds to impurities 
that locally attract or repel the interface (as for the wet- 
ting problem). In contrast, random-field (RF) disorder 
describes pinning energies which are affected by the im- 
purities inside the entire region delimited by the interface 
(Fig. [T]). This situation is encountered for example in 
magnetic systems when impurities modify the local mag- 
netic field. The disorder potential Edis generated by the 
impurities in these two cases is given by 



{V{h{x),x) random bond 

j^^^"^ duV{u, x) random field. 



(2) 



where for each value of a;, the function V{u) is a short- 
range-correlated Gaussian noise: 



V{x, u) = 0, V{x, u)V{x', u') = 6'^{x-x')R{u~u') (3) 

where the over-bar stays for the average over all disorder 
realizations and R{u) is rapidly decaying function. 

A third interaction appears in the presence of an ex- 
ternal force. With strong intrinsic dissipation, the fol- 
lowing over-damped zero-temperature equation of mo- 
tion applies: 



r]dth 



5h 



+ f = cdlh + Fd,,{h,x) + f, (4) 



where 77 is a friction coefficient and i^dis is the pinning 
force. For small external force /, the manifold ends up 
pinned in a metastable state whereas it moves with finite 
velocity v at larger force. The moving phase is separated 
from the pinned regime by the critical force /c [25[ . 



B. Self-afHne reference states (T = 0) 

Three reference states are present in the velocity-force 
diagram at zero temperature (Fig. [J): the equilibrium 



1. Equilibrium (f = 0) 

The large-scale properties of elastic manifolds at equi- 
librium are independent of temperature (26l . [27l . [3^ . [33l . 
d^l, because the disorder remains relevant. Indeed, 
sample-to-samplc fiuctuations of the ground-state energy 
grow with the system size L as Egg ~ with a positive 
exponent 9. 

The values of ^oq and 9 depend on the dimension d, 
the range of the elastic interactions and on the type of 
disorder (RB or RF), but they are not independent: In a 
system of size L, the displacement h scales as L''"^, and 
the short-range elastic energy as L^'''=<i+'*~^ (see Eq. ([1])). 
At equilibrium, the elastic energy and the disorder con- 
tribution should scale in the same way. This implies the 
scaling relation 



= 2C 



cq 



d-2. 



(5) 



The motion is governed by the minimal energy barriers 
between metastable configurations. The barrier between 
configurations which differ over a size I is of order U ^ l'^ 
for large Z, with a positive barrier exponent tp. This power 
law is responsible for the logarithmically slow relaxation 
towards equilibrium 35, 36]. The relation ip = 9 is widely 
accepted [S^l (see, however, [s^, where ijj — d/2 is pro- 
posed) . 



2. Depinning (f — fc) 

At the depinning threshold / = f^, the moving mani- 
fold is sclf-affine both in space and in time. This means 
that displacements h{t) — h{0) on scales x and at time 
intervals ~ t are statistically equivalent to displacements 
on scales aCdop/zdop^. times ~ at (here ^dop is the dy- 
namic exponent). We note that in equilibrium the loga- 
rithmically slow dynamics prevents self-affinity in time. 

Above the depinning threshold (/ > fc), the velocity 
vanishes with a characteristic exponent, v ^ {f ^ fc)^ , 
and the motion is characterized by avalanches of a diver- 
gent typical size ^, 



(6) 



In the language of critical phenomena, the velocity plays 
the role of the order parameter and the force the role of 
the control parameter (l^ . 
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The exponents (3, Cripn ; Vc \m, and Zdop are constrained 
by scaling relations [33, EO] ■ The velocity of the mani- 
fold is related to the characteristic time of an avalanche, 
t ~ ^^dop^ a^jjj ^}jg distance the manifold advances 
during this time ^^^op^ as -y ~ ^Cdop-^dop^ This yields a 
hyperscaling relation 



P ~ ^'dcpC^dop ~ Cdcp)- 



(7) 



Depinning exponents depend on the dimension d of the 
manifold and on the range of the elastic interactions 
but they are independent of the type of disorder (RF 
or RB), mer ging the two equilibrium universality classes 
into one f4ll. l42|. 

Another scaling relation holds, both for depinning and 
equilibrium, if the equation of motion preserves the sta- 
tistical tilt symmetry (STS) [IJ. The tilt is a static force 
e(x) with vanishing spatial average. If one adds a tilt 
to Eq. dH), and changes variables as h{x,t) h{x,t) + 
V~^e(x) or, in Fourier space, h{q,t) — v h{q,t) +q~'^e{q), 
the same equation of motion is recovered with a new re- 
alization of the same disorder. The statistical properties 
of the manifold are unchanged an d the response function 
behaves like in the pure system dh{q,t)/de{q) ~ q~^. 
Dimensional analysis of Eq. ^ implies that the force e 
scales with the exponent — l/i^. On the other hand, self- 
afHne displacements h scale with the exponent ^, so that: 



(8) 



2-C 



Statistical tilt symmetry symmetry is violated by a non- 
harmonic elastic energy or by certain anisotropics of the 
random medium. This violation is relevant at depinning 
where a new universal behavior is observed [2^, 0, |3l ■ 
A phenomenological mapping has been proposed between 
a one-dimensional str ing at the depinning transition and 
directed percolation [i^, |4^ . A concrete model belong- 
ing to this class contains clastic energies stronger than 
harmonic, for example 



Eei[h] 



(9) 



with C4 > [43[ . We note that the violation of statistical 
tilt symmetry changes the universality class at depinning, 
but not in equilibrium. 



3. Fast flow (f > fc) 

In the fast-flow reference state, for f ^ fc, the 
quenched pinning force reduces to an annealed stochas- 
tic noise [40] because in the co-moving frame, one has 
Fdis{h,x) = FdisiSh + vt,x) ~ Fdisivt,x). For short- 
range correlated pinning force, the strength of the disor- 
der plays the role of an effective temperature Toff, since 



FdUvt,x)Fdis{vt',x') ~ A{0)S{vit - t'))S{x - x') 

= ^sit-t')Six-x'), 
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FIG. 3: (a) Zero-temperature dynamical phase diagram, (b) 
Schematic structure factor for f > fc- The correlation length 
^ fixes the crossover between depinning and fast flow. 



where A(0) measures the disorder strength. Since 
dv/df ~ r]~^ in this regime, it follows that Teff ~ 
A(0) / (rjv) from a generalized fluctuation-dissipation the- 
orem. For the string, the fast-flow state corresponds to 
the random walk. 

Critical exponents for the three self-affine reference 
states are summarized in Table |I| 



III. DYNAMIC PHASE DIAGRAM 

At non-zero temperature, the manifold moves with fi- 
nite velocity for all forces / > 0, and the long-time dy- 
namics reaches a steady state. Possible steady states are 
contained in a dynamical phase diagram exhibiting, at 
different length scales, the three reference states of Sec- 
tionini Above threshold, this diagram is well understood 
in terms of the analogy with second-order phase transi- 
tions. 

The steady-state dynamics below threshold can be 
studied by means of our powerful algorithm which ex- 
ploits the special properties of the dynamics (see Sec- 
tion lIVp . This algorithm yields a phase diagram which 
differs from the standard picture of a second-order phase 
transition. 



A. Above threshold (/ > fc) 

We first consider the region above the critical force 
where a steady state exists even at zero temperature. 
Here the analogy of the depinning transition with crit- 
ical phenomena is well understood 



14l |: As in a ther- 



modynamic second-order phase transition, the connected 
two-point correlation function of the order parameter is 
characterized by a correlation length which diverges at 
the critical point. The steady-state velocity is the order 
parameter of the depinning transition and its two-point 
correlation function, 



{{v{x,t)~v){v{0,t)-v)) ^e-^'-'^/^, 



(10) 



indeed diverges for f ^ f^ ■ Here, brackets stand for the 
thermodynamic (or steady-state) average. In practice. 
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FIG. 4: (a) Dynamical phase diagram for / < fc for vanish- 
ing temperature, (b) Schematic structure factor for f < fc- 
The length Lopt fixes the crossover between equilibrium and 
depinning. 



this correlation function is not easily accessible because 
the steady-state velocities necessitate the long-time inte- 
gration of the equation of motion (47| . 

The correlation length ^ separates two length scales in 
the manifold (sec Fig. [3^): on scales smaller than ^, the 
geometry of the interface is characterized by the expo- 
nents of the depinning reference state (the critical phase 
in the language of magnetic transitions). In contrast, on 
length scales larger than ^, the interface is governed by 
the exponents of the fast-flow reference state (analogous 
to the ferromagnetic ordered phase of a magnetic tran- 
sition). The length ^ can be measured [47[ through the 
structure factor [48j . 



Siq) = 



1 



^ ,,„ , d'^x h{x,t) e-'9^ 
d'^x e-'«^ {h{x,t)h{Q,t)), 



(11) 



where the second equality makes use of spatial translation 
invariancc. For inverse lengths q belonging to a self-affine 
regime with a single roughness exponent C, the structure 
factor takes the form 

5(<z)~g-('^+2C). 

The crossover between the depinning and the fast-flow 
regimes can be conveniently extracted from the change 
of slope of the structure factor S{q) (see Fig. IJb). 



B. Below the depinning threshold 

Below the depinning threshold (/ < /c), at zero tem- 
perature, the manifold is permanently pinned. We first 
discuss the phase diagram in the limit of vanishing tem- 
perature, obtained with the methods of Section [TVl 

As in the regime above threshold, the structure fac- 
tor allows us to access the self-affine regimes present in 
the interface. Our results from Section |V] are summa- 
rized in Fig. a crossover length, Lopt, associated to 
the maximal barrier encountered across the optimal path 
along the system, separates two roughness regimes: On 
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FIG. 5: (a) Dynamical phase diagram at finite temperature, 
(b) Schematic structure factor below fc. The correlation 
length 5 fixes the crossover between depinning and fast fiow, 
and the length Lopt the crossover between equilibrium and 
depinning. 



length scales smaller than Lopt the roughness of the in- 
terface is described by the equilibrium exponent , cor- 
responding to the paramagnetic phase in the language 
of magnetic transitions. For distances bigger than Lopt, 
the roughness is described by depinning exponents (that 
is, the exponent of the critical phase). This is at vari- 
ance with standard critical phenomenon, where the criti- 
cal phase appears at large length scales only at the critical 
point. Our results are summarized in the phase diagram 
of Fig. [3^. Increasing the external force, Lopt decreases 
and when the depinning threshold is reached (/ f^), 
Lopt coincides with the Larkin length (sH . 

Our algorithm cannot access very small forces, but 
our results are compatible with predictions in the creep 
reg ime, that is, at low temperatures for f ^ fc [13, El, 
I49L [sol, [Slj. In this regime, a phenomenological scaling 
argument suggests that the mean velocity is produced by 
activated jumps on a length scale which diverges in the 



equilibrium limit / ^ as L, 



opt 



f- 



the typical energy barriers scale as 

uif) ~ Lopt(./)''' ~ /-^ 



At this length. 



(12) 



where = tpi"^^ is the creep exponent. The velocity-force 
characteristics for small / is thus a stretched exponential. 
The functional renormalization group [4^ predicts that 
on scales below Lopt the system is in equilibrium and that 
scales larger than Lopt are characterized by deterministic 
forward motion. 

The structure factor at finite temperature is sketched 
inFig.lSJ). It presents three roughness regimes: for length 
scales smaller than Lopt, the interface is in equilibrium. 
Between Lopt and ^, it is characterized by depinning 
exponents, and for larger lengths it is in the fast-flow 
regime. At finite temperatures, the two crossover lengths 
Lopt and ^ » Lopt diverge as / ^ (see Fig. [5^). While 
^(/, T) diverges for T — > below the depinning thresh- 
old (see Fig. [3^), we find that Lopt(/, L) saturates to a 
finite value Lopt(/, L = 0). The phase diagram of Fig. [5^ 
exhibits the connection between depinning and equilib- 
rium. 



6 




FIG. 6: Sequence of metastable configurations FIG. 8: Ground-state, critical and dominant configurations 

{. . . , Qfe, Ofe+i , . . . } detected by our algorithm. For clarity, below fc. An excitation of size Lopt escapes the system from 
afc+i is advanced by a small amount with respect to au- the dominant configuration a. The configuration /3 relaxes to 

the metastable state 7 via an avalanche of size Lroiax. 




path 



FIG. 7: Escape path from a metastable configuration via 
a configuration /3 that relaxes towards a new metastable con- 
figuration afc+i of lower energy. 



IV. T ^ DYNAMICS FOR FINITE SAMPLES 



In this section, we discuss the detailed properties of the 
dynamics below threshold in the the zero-temperature 
limit for a finite sample. We call this the 'Arrhenius 
limit', because the time At to overcome an energy barrier 
U is governed by the Arrhenius formula At ~ exp [U /T]. 
We show in the present section that in this limit the 
steady-state dynamics is characterized by a forward- 
moving sequence of metastable states of decreasing en- 
ergy, which we are able to compute using the exact algo- 
rithm described in Appendix [XI Moreover, below thresh- 
old, for each value of the external force, a unique dom- 
inant configuration is occupied with probability one, in 
the same way as for a finite system at equilibrium the 
occupation probability is entirely concentrated on the 
ground state. 



A. Ordered sequence of metastable states 

We first consider two metastable configurations a and 
7 at a fixed force /. A path Pa-ty connecting a with 
7 is a sequence of configurations labeled with the time 
parameter t, as prescribed by the dynamics of the system, 
for example a discrete rule on a lattice, or a continuum 
Langevin prescription. The path barrier Pa^-y between 
a and 7 is given by 

B[Pa^^] = m&xE[a{t)] - E\a\, (13) 

and the barrier between the two metastable configura- 
tions, -Bq^^, is defined as the minimal path barrier over 
all paths connecting a with 7. In the Arrhenius limit, 
the escape time from a metastable state a is dominated 
by the minimal barrier [E'^^'^'\a^ connecting a to a con- 
figuration with lower energy. 

^''"'^[a] = minSc-.-., with £[7] < E\q\. 

1 

The previous properties suggest the definition of a 
coarse-grained dynamics characterized by a sequence of 
metastable states, "1, ctfe of decreasing energy con- 
nected by the minimal barriers (see Fig. [5] and 
Fig. [7|). The transition from au to OLk+\ is irreversible. 
This does not mean that from a^+i, the microscopic 
dynamics cannot visit again, but rather that es- 
cape paths from a^+i towards a new configuration afc+2 
with yet lower energy exist. In fact, as shown in Ap- 
pendix 1^ the escape path can be implemented entirely 
with forward- moving steps, in the same way as for the 
zcro-tcmperature motion relevant in the dcpinning prob- 
lem [5^ Issj . The same has been pointed out |54 | for 
the Sinai model, where the effective dynamics involves 
forward moves only. 

The sequence of metastable states has two remark- 
able properties. First, all configurations situated between 
configurations and a^+i have smaller escape energies 
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than the configuration itself. Second, if in the back- 
ward direction, no metastable configuration exists which 
lowers the energy of a, then the coarse-grained dynam- 
ics starting from a is always forward-moving. These two 
properties (proved in Appendix lA 21 and Appendix lA 3p 
imply that the steady-state coarse-grained dynamics is 
always forward-directed. Furthermore, the metastable 
configuration characterized by the largest Ease belongs 
to the sequence of configurations of the coarse-grained 
dynamics. This is the dominant configuration which, as 
discussed above, is occupied with probability one in the 
Arrhenius limit. Finally, the steady state of a sample 
with periodic boundary conditions describes a periodic 
trajectory of metastable configurations, which is inde- 
pendent of the initial configuration. These properties are 
analogous of the one-dimensional problem of a particle on 
a ring, which has been solved exactly [s^, [5^ . 

The details of the algorithm arc given in Appendix [X] 
It enumerates a complete set of dynamically relevant con- 
figurations. For driven manifolds, this approach is sim- 
pler than in general [l^ because the excited configura- 
tions differ from the metastable configuration on a length 
scale Lopt which remains finite for L — > cx). In addition, 
we can restrict ourselves to forward moves. 

In Fig. [5] we compare the characteristic configurations 
of the low temperature dynamics: the dominant configu- 
ration, the ground-state, and the depinning critical con- 
figuration. These three configurations are the dominant 
states for / = 0, < / < /c and / = /c, respectively. 
Our algorithm accesses also the configurations /3 and 7 
sketched in Fig. [5] configurations a and (3 differ on a 
length Lopt, and a and 7 differ on a length irciax- ^opt 
and Lreiax are the characteristic dynamical lengths in the 
Arrhenius limit. 



V. NUMERICAL RESULTS 

Our algorithm determines, for each sample, the dom- 
inant configuration, the escape barrier i?osc, the size of 
the optimal thermal excitation Lopt, as well as the size 
of the deterministic avalanche Lreiax (see Fig. [5]). These 
results are also compared with the results obtained by 
other methods at finite and at zero temperature. 



A. Analysis of the dominant configuration 

We first consider the random-bond elastic string with 
statistical tilt symmetry. Complementing results of [2]| . 
we show in Fig. [9] the collapse of the structure factor 
of the dominant configuration at different forces, with 
all lengths rescaled by the disorder-averaged size of the 
thermal excitation, Lopt , which is obtained directly from 
the simulation. The remarkable quality of the collapse 
(which is free of adjustable parameters) is also due to the 
fact that the control parameter ///c, for each sample, is 
obtained with the sample-dependent critical force. 
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FIG. 9: Rescaled structure factor of the harmonic string in 
the Arrhenius limit (averaged over 1000 disorder realizations) 
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FIG. 10: Steady-state structure factor for a string with hard 
metric constraint in the Arrhenius limit (averaged over 10000 
disorder realizations) for L — M = 64. Curves for different 
forces are shifted for clarity. 



In Fig. [9l the change of regimes between the small-g 
region (governed by depinning) and the large q regime 
(dominated by equilibrium) is manifest. We note that 
Lopt characterizes the steady-state dynamics below fc in 
the Arrhenius limit in the same way as ^ does above the 
depinning threshold [ij, ■ 

We now consider the model with a hard metric con- 
straint, which violates STS. This case is particularly in- 
teresting because the equilibrium roughness exponents 
are identical, whereas the depinning exponents differ 
strongly. Moreover, in one dimension, the STS depin- 
ning exponent is unphysical (Cdcp ~ 1-25 > 1), and will 
probably not be observable in nature. In Fig. 1101 we 
show the structure factor for the model with hard metric 
constraints as a function of the force, from the statics, 
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FIG. 11: Steady-state structure factor for a tilted string with FIG. 12: Mean escape barrier U from the dominant configu- 

hard metric constraint in the Arrhenius limit (averaged over ration for samples L x L^"^"? with L = 32, 64, 128 (the line is 

1000 disorder realizations) for L = M = 512, 256. Curves for a guide to the eye), 
different forces are shifted for clarity. 



/ = 0, to the depinning, / = fc- The small-q region is 
again consistent with the depinning roughness exponent, 
and at large-g, we recover once more the equilibrium be- 
havior. However, the collapse of the curves for different 
values of f /fc is not possible because the two exponents 
are now very close in value. The crossover between the 
two regimes no longer depends on the length scale Lopt 
alone, but also on the microscopic parameters for which 
the string senses the presence of a nharmonic co rrections 
to the elastic energy, that is, when (/li+i — hi)^ ~ 1. 

In the absence of STS, the roughness exponents at 
equilibrium and at depinning are very close (Cdcp 0.63 
against Ccq = 2/3). However, the physics is very differ- 
ent in the two regimes and reliable signals of the no-STS 
depinning regime are present when a tilt is applied. A 
tilt can be realized through a shift s ■ L imposed on the 
interface boundary conditions (i.e. ~ ho + s ■ L and 
< s < 1). Following Ref 'ST?! the roughness exponent for 
a tilted interface will differ from that of an untilted one, 
and for a tilted string it is expected to be Ctiit = 0.5 [ssj . 
In Fig. [TT] we show the structure factor for the model 
with hard metric constraints in presence of a tilt s = 0.5. 
In the statics, / = 0, the tilt has no effect. For / > 
wc observe that tilted interfaces, at large length scales 
become less rough than untilted ones, in good agreement 
with Ctiit « 0.5. This analysis confirms that in the creep 
regime the large scale structure of the string is described 
by deterministic processes belonging to the correspond- 
ing depinning universality class of the system, regardless 
of the violation of the STS symmetry. 



B. Barriers 

In the Arrhenius limit, the velocity of the interface 
depends only on the barrier height. Assuming a narrow 
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FIG. 13: Distribution of U for 11000 disorder realizations 
(samples L x L'^^op ^^^y^ i ^ 8,16,32). The width of the 
distribution decreases with increasing L. 

distribution of barriers, we can relate the mean velocity 
of the interface to the typical barrier U{f): 

vif)=Lill(^^Y^\-^^(f^ (14) 

In the creep regime, when / is very small, we have 
ireiax ~ ^opt and C/(/) is given by Eq. dH]). The phe- 
nomenological expression for the velocity at low forces is 
the stretched exponential of the creep formula [12, El, 

v{f)^cM~P{fclfY]- (15) 

Our algorithm yields the escape barrier from the domi- 
nant configuration. For a proper scaling (59j of the sam- 
ple dimensions L x L^'*°p , we may identify this barrier 
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FIG. 14: Finite-temperature Steady-state structure factor 
from Langevin simulations of an harmonic string of size 
(L,M) = (1024,2048). The inset shows S{q) for T = 0.05 
and ///c = 0.8,0.85, . . . 1.2. The crossover length ^ is a fit 
parameter for the collapse in the main panel. 

with the characteristic barrier U{f) and test the phe- 
nomenological argument of Eq. (fT5|) . In Fig. [12] we see 
that U{f) increases with decreasing /. This is consis- 
tent with the phenomenological arguments leading to 
Eq. p2)) . even if we cannot determine the exponent fi, be- 
cause of the shortcomings of our algorithm at very smaU 
driving forces. 

In Fig. I13i we show the barrier distribution as a func- 
tion of the size L for samples scaled properly as i x i^'^'^p . 
The decay of the distribution for large U appears faster 
than exponential. Clearly, the distribution becomes nar- 
rower as the system size is increased. Extrapolation of 
these results for an infinite system leads to a well-defined 
value of U{f), confirming the phenomenological assump- 
tion of a typical barrier size. 



FIG. 15: ^ as a function of / for different T 
increasing from the top to the bottom (T — 
0.0025, 0.005, 0.01, 0.025, 0.05, 0.075, 0.1, 0.15, 0.2). 

all data at small q by using the scaling S{q) / S{2ti / ~ 
with a fitting parameter ^. 

The collapse of Fig. [14] provides the correlation length 
^ as a function of force and temperature, and wc show 
it in Fig. [15] as a function of the reduced force ///c for 
different T. ^ decreases monotonically with /, and tends 
to diverge near in the T — ^ limit, showing that ^ is 
ultimately controlled by the velocity of the string. This 
behavior is fully consistent with the schematic phase di- 
agram of Fig. [5] 

In summary, the Langevin simulations confirm the sce- 
nario valid for low-tcmpcraturc Arrhcnius dynamics for 
finite T . We find no indication of a divergent (or even 
increasing) correlation length as we approach /". The 
main modification at finite T is that ^ diverges at / = 
instead of at /c, since this divergence is controlled by the 
vanishing of the steady-state velocity. 



C. Finite temperature 

Our algorithm of Section IIVI is powerful, but it takes 
the T limit before the thermodynamic limit. For 
a macroscopic physical system, the order of these lim- 
its should be interchanged. To confirm that the order of 
limits does not affect physical properties, we have thus 
performed Langevin simulations of the equation of mo- 
tion Eq. ^ at small finite temperatures with parameters 
•q = c = 1 and with random-bond disorder from a normal 
distribution interpolated with a cubic spline, as in [60j . 

In the inset of Fig. [T31 we show S{q) for the steady- 
state motion of the string at T = 0.05 for different forces 
around fc- Two regimes are always present. They cor- 
respond to the depinning at small scales and the flow 
regime at large scales. The crossover between these 
regimes corresponds to the correlation length ^(/) of 
Fig. [5] In the main panel of Fig. [TJ] we have collapsed 



D. Deterministic avalanches below fc 

In the previous sections we considered two length scales 
below threshold, namely the size Lopt of the optimal ther- 
mal excitation, and the correlation length ^ above which 
the disorder acts effectively as a thermal-like noise. Both 
lengths appear in steady-state quantities, such as the 
structure factor S{q) or the velocity v of the interface. 

Another length scale, irciax can be identified as the 
typical size of the deterministic avalanche that drives the 
interface from the optimal activated jump of size Lopt to 
the next metastable state. 

Lroiax thus mcasurcs the distance between consecutive 
metastable states. As shown in Fig. [TO] Lrdax ~ Lopt 
diverges approaching /c from below, with the character- 
istic exponent fdop. However, Lrciax does not describe 
steady-state properties below the depinning threshold, 
and represents no genuine divergent length scale below 
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FIG. 16: Behavior of the dynamical length Lreiax- This length 
diverges with the exponent ^dcp (see inset). 

the dynamic phase transition. 

A length scale analogous to Lreiax appears also in the 
transient dynamics at T = for an initially flat configu- 
ration which relaxes up to the first pinned metastable 
state [H, |6l|. Here, a crossover separates the large 
length scales keeping memory of the initial condition 
from the short length scales characterized by the dcpin- 
ning roughness exponent. This length scale diverges with 
the same exponent Vdcp as £roiax- Wc can identify it with 

-^rclax -^opt ■ 

VI. DISCUSSION AND CONCLUSIONS 

We have studied in this paper the low-temperature dy- 
namics of a driven elastic string in a disordered medium 
below the depinning threshold. In the (Arrhenius) limit 
of vanishing temperature, and at large scales, the string 
behaves as at the critical force (/ = /c), and it shows the 
universal properties of the depinning transition. This 
result contradicts the quasi-equilibrium picture of the 
creep motion, which assumes that a small force simply 
moves the system from one metastable equilibrium state 
to another, without changing the geometrical properties 
of these states [i^ . It also demonstrates that the analogy 
of the depinning transition with equilibrium critical phe- 
nomena is incomplete. The effect of even a small force 
(/ ^ fc) is in fact more radical: it drives the system away 
from the equilibrium and into the depinning regime. The 
equilibrium behavior of the string is restricted to small 
scales, The crossover between the depinning and the equi- 
librium regimes takes place at the length scale iopt of 
the optimal barrier. At the critical force, Lopt does not 
diverge, and it equals the corresponding Larkin length. 
Lowering the force, Lopt increases, diverging in the limit 
of vanishing force. Our findings are compatible with the 
FRG predictions 

To arrive at these conclusions, we have used a pow- 



erful algorithm, which is exact under the condition that 
the T —> limit can be exchanged with the thermody- 
namic limit. We have also checked its conclusions with 
conventional Langcvin simulations. At finite tempera- 
ture, a much larger length scale ^ is introduced by the 
finite velocity of the system. Length scales above ^ are 
described in terms of the universal fast-flow regime. 

The approach of the critical force is not characterized 
by a divergent length-scale of the steady-state proper- 
tics. Nevertheless, a divergent length-scale Lrciax at fc 
describes the typical size of the deterministic avalanches 
that are triggered by the activated events. This deter- 
ministic part of the motion does not affect the steady- 
state geometry and can for example not be identified in 
snapshot of the line's positions. 

The distribution of barriers is narrow and the typi- 
cal barrier grows with decreasing the drive, in agreement 
with phenomenological arguments. This shows that the 
motion is controlled by the typical barriers, yielding a 
finite steady-state velocity, rather than by rare barriers 
described by extreme value statistics j63 |. 

We conjecture that our conclusions remain valid for d- 
dimensional manifolds moving in d + 1-dimensional space 
with short-range or long-range elasticity. Moreover, the 
violation of the statistical tilt symmetry changes only 
the large-scale geometry according to the change of the 
universality class at the depinning transition. This may 
be relevant for the interpretation of creep experiments in 
thin magnetic films [l|,13|, since systems violating STS 
display a depinning roughness which is very close to the 
equilibrium one. 

Our results may apply to different experiments probing 
the creep motion of elastic interfaces. We expect that an 
experimental verification of our results is possible using, 
for instance, imaging techniques for magnetic [11, |3, \^ or 
electric d, @ domain walls in thin films: iopt could be 
extracted from the analysis of a spatial correlation func- 
tion, and Lrciax could be measured by comparing consec- 
utive (long-lived) metastable states when / is close to fc, 
since then Lroiax controls the distances between succes- 
sive metastable states. Due to its transient nature Lrciax 
could be also measured by transient methods at T = 0, 
where wc relax an uncorrelated initial condition until it 
locks to the first metastable state at a given force. 

Steady-state noise measurements (such as acoustic 
emission noise) below the threshold could also provide 
an indirect verification of our results. We expect the 
force dependent barriers shown in Fig. [T^ to control the 
waiting times between events associated with determinis- 
tic avalanches of a diverging size Lroiax, translating into 
large noise peaks. This situation could be realized in 
fcrromagnets [6^ or in material failure [gJ, [b^ . 

We thank P. Le Doussal and K. J. Wiese for illuminat- 
ing discussions all along this work. This work was sup- 
ported in part by the Swiss NSF under MaNEP and Di- 
vision 11. 
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APPENDIX A: ALGORITHMIC DETAILS 

Our algorithm computes the dominant configuration 
of an elastic string moving on a two-dimensional discrete 
L X M lattice with periodic boundary conditions both in 
L and in M. The line is described by the variables h{i), 
giving the displacement of the string on the slice i, with 
< i < L. The energy of the line is given by 



f/f=0.5 



f/fc = 0^6 



E ■ 



E ^(M* + 1) - K^? + V{^, h{^). (Al) 



Periodic boundary conditions in M are accounted for 
by the periodicity of the disorder potential: V{i, h) = 
V{i, h + M). Besides the harmonic elastic energy (which 
preserves STS), we also consider a hard metric constraint: 



\h{i) - h{i - 1)1 < 1 



(A2) 



which violates STS. 

We use elementary moves of the 'variant Monte Carlo' 
(VMC) algorithm of Ref. gl, [s^ which allows for the 
simultaneous motion of A:+ 1 adjacent sites by one lattice 
spacing if no move of k sites is energetically favorable. 
This choice of dynamics avoids certain pathologies of the 
single-site dynamics [H, [52| . 

For each sample of the disorder potential, the equilib- 
rium ground state [1^, [sl] , the critical depinning force /c 
and the associated zero-temperature configuration at 
[isl . [5^ can be computed easily. 

Below threshold, two kinds of motions are present: 
the first is deterministic, as defined by the VMC algo- 
rithm, and it relaxes each unstable configuration toward 
a metastable state. The second is the activated dynam- 
ics connecting a^. to afe+i through the sequence of VMC 
moves belonging to the optimal path, which is character- 
ized by the barrier Ecsc[oik\- 



1. Complete-enumeration scheme 

We initialize the dynamics from the equilibrium 
ground-state, which by its very nature cannot relax 
through backward moves when the external force is pos- 
itive. This allows us to restrict our attention to forward 
moves only (see Appcndix IA 3|) . At a given external force, 
we let the ground state relax toward the metastable con- 
figuration ofQ. For each transition au ^ otk+i we build 
the archive of the visited configurations, [3i, with increas- 
ing energy E[(3i\ < E[(32] < E[(3:i] < A constant E^^t 
is also introduced in order to compute Eesc\ctk]- The 
archive is initialized with a single configuration, f3i = ak, 
and -Ecut is set to zero. At each step the configuration 
/3i is taken and erased from the archive. Two operations 
are performed on this configuration. First, we update the 
maximal barrier 

Ecnt = max(£'cut, - E[ak]) 




f=o. 



Jo 



///e = < 




FIG. 17: Metastable states of the low-temperature in a 
{L,M = 32,64) system, for different forces / < /c. The 
dominant configuration is emphasized. 



Second, wc relax /3i under the VMC dynamics to a con- 
figuration /3'. All configurations connected to (3' by VMC 
moves are incorporated to the archive. If has lower en- 
ergy than ak , the construction ends and the barrier Ecut 
as well as the last configuration erased from the archive 
f3i and its associated metastable state, are output. 

We can identify all the quantities defined in the previ- 
ous sections: 



ak+i = /?' 

Ecsc[0ik] — Ecut 



(A3) 



where Q{x) is the step function with the prescription 
0(0) = 0. Let us remark once again that for each real- 
ization of disorder we analyze the data corresponding to 
the configuration with the maximum value of Ease- In 
Fig. [T7]wc sketch the output of our algorithm for a small 
sample at different driving forces. 



2. Bound on escape energies 

In this section, we show that all configurations be- 
tween the metastable configurations ak and ak+i have 
smaller escape energies than the configuration ak itself. 
This assures that the coarse-grained dynamics must pass 
through the metastable configurations with the largest 
barrier. 
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FIG. 18: Decomposition of a configuration a{t) as a sum of 
ai{t) and 02 (t), the configuration (f> by the sum of (j>i and <^2, 
the configuration a^(t) by the sum of 01 and a2(t), and the 
configuration <j)a{t) by the sum of ai(t) and (fi2- 



Let US consider a configuration a which moves to 7 in 
order to pass through the minimal escape barrier E^'^. 
We show that any configuration with h^[x) > ha{x) 
(Vx) and with an energy barrier E'^^ > E^'^ satisfies 
hci){x) > /i-y(x) (Vx). 

We suppose that the property is not true which means 
that hfi,{x) < hj{x) for some x. We follow the dynamical 
evolution of the line a(t) and compare it with the evolu- 
tion of two sets of configurations. The first one, a^{t), is 
defined as 

h.jMx) = h^''^ '^h.it){x)>h,{x) ^ ^^^^ 
\ha(t){x) otherwise 

The second one, <j>a{t), is defined as 

, f,)h4,{x) a ha^t){x) < h^{x) 
\'^a{t)[x) otherwise 

We start by comparing a{t) and <f)a{t)- Because of the 
metastability of (j) we expect that, initially, E[a^{t)] < 
E[a{t)]. However, this situation cannot continue up to 
tend because, by definition, a(t) is the path which en- 
counters the smallest barrier. We call i the smallest time 
at which 



E[a{i)] < E[a^{i)]. 



(A6) 



Using the decomposition of Fig. [181 for t = i, we can 
write 

E[a{i)] - E[a4,{i)] = 

E[ai{i)] - E[Mi)] + ^E,i[a^{i) ^ a{i)], 

where the variation of elastic energy is understood as 

AE,i[a4,{i) ^ a(i)] = E,i[ai{i),a2ii)]~ E,i[Mi),a2ii)]. 



Let us now compare the energy associated to the config- 
uration (t) to the energy associated to the metastable 
state (j): 

E[^^{i)]~E[^] = E[a,{i)]~E[Mi)] + AE,,[(f> ^ 

The convexity of the elastic energy assures that 

AE,i[<j) ^ < AE.ila^i) a{i)]. (A7) 

Using Eq. (jA6| we conclude that E[(l)a{t)] < £^[0]. All 
paths connecting with (pait) niust overcome at least 
a barrier E^^. This means that there exists a time t* 

{0 <t* < i) for which -E;[(/>a(i*)] - E[(j3] > E"^''. We can 
use the usual energy decomposition also for the configu- 
ration visited at the time t* (see Fig. [iH] for t^t*). The 
following inequality it is easy to prove: 

E[dp^{t*)] ~ E[4>\ < E[ait*)] - E[a^{t*)]. (A8) 

This implies 

E[a{t*)]>E[a^{t*)]+El^^. (A9) 
On the other hand 

E[a] < E[a^{t*)]. (AlO) 

It follows that 

E[a{n]-E[a]>Ef'^-. (All) 

This contradicts the initial assumption, demonstrating 
that all configurations between the metastable configu- 
rations Q!fc and Qffc+i have smaller escape energies than 
the configuration ak itself. 

3. Effective forward motion 

In the present section we prove the following state- 
ment: If there is no metastable configuration which low- 
ers the energy of a in the backward direction, the coarse- 
grained dynamics starting from a will be forever forward- 
directed. In other words, there exists an escape path 
starting at a involving only forward motion. This al- 
lows us to restrict the search of new configurations, and 
reduces the complexity of our algorithm. It also allows 
us to discard the archive of accumulated configurations 
whenever a new metastable minimum configuration is en- 
countered. 

We suppose that the above statement is false by as- 
suming that the configuration 7 relaxes towards a config- 
uration (j) such that, for some x, we have h^(x) < h^{x). 
Due to the convexity of the elastic energy we may restrict 
ourselves to the region where the relaxation is strictly 
backward and h^{x) < h^{x) for all x. The stability of 
a with respect to backward movements imposes 



ha{x) < h^(x) < hj{x) \/x. 



(A12) 
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The energies associated to these eonfigurations satisiy 
E[a] > E[j] > E[(l)]. (A13) 



We now show that given Eq. (|AT2l) and Eq. ((Al3l) . there 
exists a configuration 0i for which 



E[a] < E[-f] > E[(f>] > E[(l)i] 



(A14) 



In analogy with Appendix lA 21 we can compare the evo- 
lution of a{t), a^{t) and (pait) defined in Eq. (|A4[) and 
in Eq. (|A5p . The configuration corresponds to 4>a{t): 
where t is the smallest time at which 



E[a{i)] < E[a^{i)]. 



(A15) 



This construction can be applied to configurations 
(/)2, 03, . . . up to 0„ = 7. In this case we have E[(/)n] = 
E[^] and the statement is shown to be correct. 

APPENDIX B: LONG-RANGE ELASTICITY 

For interfaces with long-range interactions, the elastic 
energy can be written in compact form in the harmonic 
approximation: 



E, 



d'^q\q\''hgh. 



9- 



(Bl) 



Here, the parameter a controls the range of the interac- 
tions. The standard short-range interaction corresponds 
to a = 2. The long-range interactions acting on the con- 
tact line of a liquid meniscus [l^ and on a propagating 
crack front [2^ yield a = 1. 

We expect the general scenario presented in this paper 
to remain valid for the long-range case. The numerical 
values of the universal exponent should depend on the 
range of the elastic interactions, parameterized by a. 



The scaling relations discussed in this paper can be 
adjusted to the case of general a. From a dimensional 
analysis of Eq. (|Bip we infer the generalization of Eq. ([5|) 



9 = 2C. 



(B2) 



The four exponents of the depinning transition are still 
constrained by the hyperscaling relation Eq. ([7]). In pres- 
ence of STS, the argument gi ven for a = 2 h olds and the 
response function is given by dh{q,t)/de{q) ^ q~°'- The 
STS scaling relation writes 



1 



a ~ C, 



(B3) 



A real-space representation of Eq. (|Bip involves frac- 
tional derivatives [6^. 

A discrete version of the force derived from Eq. (jBip 
is given by [i^, E^l 



(B4) 



It presents strong finite-size effects for a ~ 2 and all 
choices a > 2 corresponds to the standard Laplacian 
force. It is more convenient [6^ to use the following 
discretized force 

= (B5) 



with 



r(|i-j|-f)r(« + i) . 



where r(a;) is the Gamma function. This discretization 
holds for all a > 0, and it is unaffected by slowing down 
for a ~ 2. 
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